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Abstract 

Distribution functions defined in accord with the quantum theory of mea- 
surement are combined with results obtained from the quantum Langevin 
equation to discuss decoherence in quantum Brownian motion. Closed form 
expressions for wave packet spreading and the attenuation of coherence of a 
pair of wave packets are obtained. The results are exact within the context 
of linear passive dissipation. It is shown that, contrary to widely accepted 
current belief, decoherence can occur at high temperature in the absence of 
dissipation. Expressions for the decoherence time with and without dissipa- 
tion are obtained that differ from those appearing in earlier discussions. 
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The problem of decoherence in quantum systems has been of considerable recent interest. 
A sampling of earlier theoretical work appears in Refs. [1-5], but for a detailed survey we 
refer to a recent book devoted wholly to the subject [6]. In an introduction to the contents 
of this book, Joos surveys the current situation and, in discussing the mechanism of deco- 
herence, states that "irreversible coupling to the environment seems to have become widely 
accepted (and even quite popular) during the last decade, not least through the various 
contributions by Woljciech Zurek and his collaboratorsln this present work we come to quite 
the opposite conclusion and in fact show that, in the high temperature case considered by 
Zurek [4] and others, decoherence occurs in a characteristic time [ see Eq. (21) below] 
that is independent of the Ohmic decay rate 7 [see Eq. (19) below] which characterizes the 
strength of the coupling to the environment. Furthermore, the formulation we use is exact 
and enables us to show explicitly that previous estimates of the decoherence time arise from 
an inconsistent application of long-time asymptotic formulas to obtain a short-time result. 
Apart from its relevance to the question of classical-quantum correspondence and the foun- 
dations of quantum mechanics, this work, especially that part dealing with entangled states, 
is clearly relevant to the host of experiments on decoherence [7], quantum teleportation [8] 
and quantum information and computation [9,10]. 

Much of the discussion of decoherence has been in terms of the simple problem of a 
particle moving in one dimension that is placed in an initial superposition state (Schrodinger 
"cat" state) corresponding to two widely separated wave packets. Decoherence is said to 
occur when the long-time interference pattern is destroyed. The key questions asked are, 
first, under what conditions does decoherence occur and, second, what is the decoherence 
time. Previous discussions of the problem have either used the Feynman- Vernon influence 
functional technique [1] or master equation techniques [2-4] and have been confined to the 
case of Ohmic dissipation. In either case it is assumed that the initial state is decoupled 
from the environment. Here we use a different formulation in terms of quantum distribution 
functions introduced some time ago by Ford and Lewis [11]. With this formulation we 
are able to obtain exact closed form expressions for the spreading of a wave packet and 
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for the attenuation of interference in the two wave packet problem, which are valid for all 
temperatures and for a very general dissipative environment. In particular, we avoid the 
assumption that the particle is initially decoupled from the environment; the particle is in 
equilibrium with the environment at the time it is put into the initial state by a measurement. 
An important conclusion is that decoherence can occur for 7 — > (absence of dissipation). A 
further feature of the Ford-Lewis formulation is that the density matrix for the entire system 
is employed, i.e., one does not trace over the environment to obtain a reduced density matrix 
as previous investigations have done. 

In the formulation of Ford and Lewis, we regard the particle as part of a larger system 
of particle coupled to a reservoir. Initially (or in the distant past) the complete system is in 
equilibrium at temperature T, described by the normalized density matrix 

e -H/kT 

P° = Tr{e-^/ fcT } ' ^ 
where H is the system Hamiltonian. To a measurement of x one associates a measuring 
function a{x\) (x 1 is a c- number) such that \\a(x — Xi)$|| 2 dx\ is the probability that if the 
system is in a normalized quantum state $ the instrument will read between x\ and Xi + dx±. 
An example is the "Gaussian slit" [12] for which 

a(x 1 ) = (2nal)-^eM-^h (2) 

where <j\ is the experimental width. It follows that if x is measured at time ti, the probability 
that the instrument will read between x\ and dx\ is W{x\,t\)dx\, where 

W(x u h) = Tr{a[x(t!) - x 1 ]p a[x(t 1 ) - x^}. (3) 

Here x(t) is the Heisenberg operator at time t, 

x(t) = e iHt/n xe- iHt/n . (4) 

In the same way, if x is measured at time t\ and again at a later time £2 the probability that 
the first instrument will read in range dx\ and the second in dx 2 is W(xi,ti, X2,t 2 )dxidx2- 
Using an obvious shorthand notation, 
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W(l, 2) = Tr{«(2)a(l)p «(l) t «(2) t }. (5) 

Here in a(j) = a[x(tj) — Xj] the index j is meant to indicate not only the time tj and 
instrumental position Xj, but also the instrumental parameters such as a width Uj. In this 
way one can go on to define higher order distribution functions. 

The distribution functions can be expressed in terms of the corresponding characteristic 
functions, 

W(l) = /_"f^(l)e-**\ 
W ^2) =r^f ^h^ 2 )e-^ +k ^. (6) 

J-oo 1lX J -oo 2tv 

Now the key formulas needed from Ford and Lewis [11] is that for quantum Brownian motion 
these characteristic functions are given by the general formulas 

t(l)=e X p{-±k*(x(t 1 )*)} 

f°° dqi feu,., k lx 

j 2 2 

£(1, 2) = exp{-- £ E «tj>(ti) + <U)x{t 3 )) kjki} 
d<7i f°° dq 2 -X „ , kj 



x exp{gi/c 2 [x(ti),a;(t 2 )]}. (7) 



Here a is the Fourier transform of the function a describing the measurement, 

/•oo 

a(l) = / dx 1 a(x 1 )e- tqiXl . (8) 



We should remark that in the derivation of these formulas it was necessary to assume that 
the commutator [x(ti),x(t2)] is a c-number. This is the case for quantum Brownian motion 
[13]. 

We first apply these formulas to obtain an expression for the spreading of a wave packet. 
That is, at time t\ a measurement is made with an associated function of the form (2) and 
then at a later time t 2 a second measurement of the same form is made (with index 1 — > 2). 
The integrals are all standard Gaussian integrals, and we obtain the results 



where (note the misprint in the Eq. (7.18) of [11]) 



a~ = aj + (x 2 , , 



2apr = -s(t 2 -h)+2 (x 1 ) . (10) 
In this last, s(t 2 — 1±) is the mean square displacement 

s(t) = ({x(t 1 )-x(t 1 + t)} 2 ). (11) 
The distribution W(l,2) is a Gaussian quadratic form, with mean square width given 

by 

/OO /'OO 
dx 1 / dx 2 {x l -x 2 ) 2 W{l,2) 

2 [rr(ti),a;(ti +t)] 2 2 

= a l 7^ 1" + CT 2- I 12 ) 



This is an exact general formula for the spreading of a Gaussian wave packet, expressed in 
terms of the mean square displacement and the nonequal-time commutator. For the special 
case when the second measurement is made with infinite precision (er 2 = 0) and for a free 
particle without dissipation and at zero temperature (s(t) = 0, [x(t 1 ),x(t 1 + t)] = iTit/m) 
this reduces to the familiar formula of elementary quantum mechanics [15]. 

In the case of an unbound (free) particle the mean square displacement (x 2 ) diverges. 
We can obtain a simple expression in this limit if we introduce the conditional probability 

W(l 2) 

P^-^h- tl ) = J^-^>. (13) 



Using the expressions (9) we find 



ex P{-2^t)} 



p( X ,t)= v, z ■ (14) 

y 2llW (t) 



Thus, the conditional probability is a normal distribution with variance w 2 (t). Here we 
should refer to the work of Hakim and Ambegaokar [14], who use path integral methods for 
the special case of a free particle interacting with an Ohmic bath to obtain an equivalent 
expression for wave packet spreading, (note the misprint in their expression (38) for the 
width) 

Next, we consider the case where the initial measurement forms two widely separated 
wave packets, which corresponds to the measurement function 

e*P{-^} + exp{-feif3i} 

[8lT<T?(l + e-^«)2]l/4 ' W 

where d is the separation of the wave packets, the width of each being o\ and x\ being the 
center of the wave packet pair. The second measurement is then made with a single slit 
instrument corresponding to a function of the form (2) (with index 1 — > 2). Again, the 
integrals are all standard Gaussian integrals, and we obtain results for W(l) and W{1,2). 
Again, there is a considerable simplification if we introduce the conditional probability (13). 
We find 

P(x,t) = 



v / 2^2(l + e- d2 / 8(T ?) 

x exp{ — — L }cos — \ 9 \ 

y L 1 2w 2 J Aia\w 2 

+ r xp{ -4^ }+ 2 exp{ -4^ } j- (16) 

This conditional probability is the sum of three contributions, corresponding to the three 
terms within the parentheses. The second and third clearly correspond to the sum of proba- 
bilities of the form (14) from two single slits, while the first term (that involving the cosine) 
is an interference term. It is of interest to study the ratio, a(t 2 — ti) of the amplitude of the 
interference term to twice the geometric mean of the other two terms, which we will refer 
to as the attenuation factor. We find 

a[t) = exp{ - \aW(t) ] (17) 
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In general, the interest is in the case where the second measurement is made with infinite 
precision, so in the following discussion we set o~ 2 equal to zero. 

For quantum Brownian motion, the mean square displacement and the commutator are 
given by the formulas 

2h r°° , hm 



'2Ji roo flu 

s(t) — — / dcolm{a(uJ + i0 + )} coth — — (1 — cosut), 
7r Jo 2kT 

2ih r°° 

\x(t 1 ),x(t 1 +t)} = / dulm{a(cu + iO + )}smcut. (18) 

7T JO 



where a is the response function. In the so-called Ohmic case, where the mean motion is 
m (x) + iwy (x) = 0, 

lm{a(u + *0+)} = 1 (19) 

Consider first the case of vanishingly small dissipation (7 — > 0). Then, setting a\ — 
and putting s = ^t 2 and [x(ti),x(ti + 1)] = i^t in (12) and (17) we see that 

d 2 

a(t) = exp{ = -} no dissipation, (20) 

8a 2 + 2\ 2 + 8^ 

where A = h/VmkT is the mean thermal de Broglie wavelength. Now the interest is always 
in the case where the wave packets are widely separated, d ^> a±. From this expression we 
see that for long time, the attenuation factor will be small (i.e., there will be decoherence) if 
the temperature is sufficiently high that the mean de Broglie wavelength is small compared 
with the spacing, d ^> A. The characteristic time for decoherence to occur will then be 

- - £• < 2i > 

where v = \J kT/m is the mean thermal velocity. This decoherence time is the time for a 
particle travelling with the mean thermal velocity to traverse the slit width multiplied by the 
ratio of the slit width to the slit spacing. Thus we see that we can have decoherence without 
dissipation (in the sense that is independent of 7, which characterizes the strength of the 
coupling to the environment). 
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Next we consider the case of Ohmic dissipation at high temperature, where by high 

temperature we mean kT 3> fry. Then, using (19) we see from the formulas (18) that 

, , 2kT . l-e~ 7 \ 
s(t) = (t , 

7 

[x(t 1 ),x(ti+t)\ = — (l-e- 7t ). (22) 
m7 

For short times (7* <C 1) these reduce to the those for the case of vanishingly small dissipa- 
tion, for which the decoherence time is given by (21). Thus if 77^ = 'yaf/vd is small (and 
this will generally be the case at high temperature) the decoherence time will be the same 
as for the case of vanishingly small dissipation, given by (21). 

These exact results are strikingly different from those obtained by previous investigators 
[2-5]. It appears that the disagreement arises from the fact that others have implicitly 
used a long time (7^ ^> 7 _1 ) approximation to obtain characteristic decay times. To see 
how this comes about, we evaluate the high temperature formulas (22) for very long times 
(7* > 1). Putting the result in (17) and setting ^ = we find a(t) ~ exp{— d 2 ^t/\ 2 }. This 
is an exponential decaying with characteristic time 7~ 1 A 2 /d 2 , which is exactly twice the 
expression for the decoherence time obtained by previous authors [3-5]. But we see that it 
is inconsistent with the assumption of long time used to obtain it. At short times, as we have 
seen, we recover the estimate (21). One possible reason why others have obtain inconsistent 
results may be due to the fact that they assumed that the system and its environment 
are initially decoupled whereas, by contrast, in the formulation we use the particle state is 
entangled with the environment (i.e., in equilibrium) at the time it is put into the initial 
state by a measurement. Putting this point in another way, previous discussions have been 
in terms of the reduced density matrix but, as pointed out by Ambegaokar [16], on such 
short time scales the time evolution does not operate on the reduced density matrix alone." 

In conclusion we have seen that the simple and general formulation of quantum measure- 
ment given in [11] provides a powerful method for discussing quantum stochastic systems. 
The formalism is in terms of quantum distribution functions and, when combined with re- 
sults obtained from the quantum Langevin equation, has enabled us to obtain exact explicit 
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expressions for wave packet spreading and the coherence attenuation factor. In discussing 
the latter we have seen that decoherence occurs at high temperature with or without dissi- 
pation. In either case the decoherence time is the same, given by (21). At zero temperature, 
decoherence occurs only in the presence of dissipation. 
Note added in proof: 

In order to understand more clearly the origin of our (no dissipation) result (20), we 
recently showed that it may be derived in a simple manner solely within the framework of 
elementary quantum mechanics and equilibrium statistical mechanics [17]. In addition, we 
have recently obtained an explicit general solution of the exact master equation [18]. When 
applied to the situation considered in previous discussions, namely a particle at temperature 
zero suddenly coupled to a bath at high temperature, we are led to an expression for the 
decoherence time differing by a factor of 6 from the conventional result. We see therefore 
that the conventional result corresponds to a particle that is "warming up" over a time of 
order 7 _1 ; our result corresponds to a particle that is initially at the same temperature as 
the bath. 
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